A density functional study of pressure induced superconductivity in P and its 

implication for spintronics. 
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The stability of high-pressure phases of P has been studied using density functional theory and 
the local density approximation. Using a linear response technique, we have calculated the phonon 
spectrum and electron-phonon interaction for bcc P and predict it to be superconducting with T c 
of 19 K. We propose that this phase might be realized in epitaxial thin films using templates such 
as V(100), Fe(100) or Cr(100) relevant to spintronics applications. 

PACS numbers: 61.50.Ks, 74.25. Jb, 74.78. Fk, 85.75.-d 



The study of properties of materials under pressure is 
a key experimental technique and recent discoveries of 
superconductivity in such simple metals as compressed 
Fe 0) and Lj_|3| are currently generating enormous in- 
terest 0, 0, 0f. Black insulating phosphorus, with an 
orthorhombic A17 structure, which is the most stable 
form at ambient conditions, exhibits a series of pres- 
sure induced phase transitions. 0, 0, IS 01 The first is 
a transformation to a metallic phase with the simple cu- 
bic structure at 10 GPa. This phase is known to su- 
perconduct with T C =10K jjj. An observation of a new 
simple bcc structure above 262 GPa has been made re- 
cently , which points out that similar or even stronger 
electron-phonon interaction and pressure induced super- 
conductivity may be also realized in bcc P. Growth of 
P containing materials by thin film techniques has been 
developed very recently [l(| to fabricate the class of wide 
gap phosphide ferromagnetic semiconductors which may 
provide an ideal opportunity to stabilize the bcc phase 
of P at ambient conditions using suitable templates such 
as the V(100), Fe(100) or Cr(100) substrates. This may 
lead to a significant breakthrough in technological spin- 
tronics applications in which combined spin and super- 
conducting degrees of freedom will provide a new level of 
functionality for microelectronic devices. 

In this work we address both the issue of the phase sta- 
bility of P as well as the superconductivity of its highly 
pressurized bcc structure. We use density functional to- 
tal energy prij and linear response [l2T | techniques which 
have proven to provide a reliable description of various 
ground state properites for a large class of materials. 
We also study the electronic structure of the Fe5/P„/Fe4 
(n=3,5) super lattices, with the lattice constant of bulk 
Fc to investigate the problem of how to stabilize the bcc 
phase under ambient pressure conditions. 

In the past, structural phase transformations in P have 
been the subject of many experimental [S 0> HI El and 
theoretical [jj, Q density functional theory (DFT) 
based studies. Orthorhombic A17 structured P trans- 
forms to the rhombohedral A7 phase at 4.5 GPa.|6J The 
puckered P layers in the ambient-pressure A17 structure 
have strong covalent intralayer bonding and weak bond- 



ing between layers that result in an anisotropic com- 
pressibility. Under pressure, the distance between the 
interlayer P atoms in the A17 structure decreases faster 
than the intraplane distance forming the semimetallic A7 
structure. This A7 phase, which is similar to the struc- 
ture of such group-V elements as As, Sb, and Bi, under- 
goes a transition to the simple cubic (sc) metallic phase 
at 10 GPa. 6] In the group-V elements, this A 7— >sc tran- 
sition is observed also in Sb at 7 GPa. 

With increasing pressure, at 137 GPa, the sc struc- 
ture transforms into a simple hexagonal (sh) one.|7] This 
sc^sh transition, which is accompanied by a large vol- 
ume reduction of 7.6 %, is the first such observation re- 
ported for an elemental system. The sc-sh transforma- 
tion probably occurs via an intermediate phase |7j, and 
results in the change in the co-ordination number from 
6 to 8. Since the sh structure has a low atomic packing 
fraction, namely 0.605, a phase with a higher packing 
fraction such as bcc, fee or hep should be expected at 
higher pressures. Indeed, with increasing pressure, the 
c/o ratio of the sh phase increases from 0.948 to 0.956, 
and, finally, above 262 GPa, a new high-pressure phase is 
observed, which has been proposed to be bcc.0 A high- 
pressure bcc structure has also been found in Sb at 28 
GPa and in Bi at 7.7 GPa.0] In contrast, the sh phase 
of Si and Ge transforms to the hep structure. ^7| 

To date, despite these data, no full theoretical descrip- 
tion of the phase stability in P under extreme high pres- 
sure has been reported. There are a few DFT calculations 
on the structural stability of P under pressure. Us- 
ing the pseudopotential method, Chang and Cohen 
reproduced the Al7^A7^sc phase transitions, with the 
crystal energy of the A7 phase being shifted by 2.3 mRy, 
but failed to find any stable closed-packed P phase at 
high pressure. In contrast, Sasaki et al. 0] inferred 
a sc-bec transition at a pressure of 135 GPa, rather far 
from the recent experimental observations. [7j This paper 
shows that ab-initio DFT calculations can successfully 
describe all the structural transitions in the phase dia- 
gram in P and, in particular, quantifies the sh-bec tran- 
sition, the highest-pressure structural transformation ob- 
served up to now for an elemental material. 
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FIG. 1: Phase diagram of P at high pressures. 

FIG. 2: The Eliashberg spectral function of bcc P. 



The use of the Debye-Griineisen theory within DFT 
calculations is the simplest way to investigate structural 
phase transitions at finite temperatures. The character- 
istic Debye temperature 8 d , calculated in terms of the 
bulk modulus, can be used to constuct the free energy 
F(T,V) as the function of Qd and V. The successful 
applicability of this model, in despite of its harmonic 
approximation for lattice vibrations, has been demon- 
strated |18| for some elemental materals under high pres- 
sure. To construct the phase diagram, the thermody- 
namical Gibbs potentials G(P, T) = F(V, T) + PV are 
calculated and compared for various structures while the 
isothermal pressure dependence P(V) is obtained by di- 
rect differentiation P = —(dF/dV) T for each structure. 
We use the full potential linearized augmented plane 
waves (WIEN) code [l^, which is one of the most ac- 
curate DFT schemes. Exchange and correlation effects 
are treated using the generalized-gradient approxima- 
tions ppj | while muffin-tin radii of 1.7 a.u. are fixed for all 
structures. As Fig. 1 of the phase diagram of P shows, 
the calculations reproduce quantitatively well all avail- 
able experimental data on the pressure induced struc- 
tural transformations in P. In particular, at just above 
the critical pressure where the sh-bec transition takes 
place, at 260 GPa the bcc phase is a real possibility and 
the bcc phase remains stable at still higher pressures. 
The free energy calculation of the bcc phase shows a mini- 
mum at 3.05 A. We predict also that the sh-bec transition 
has a positive dT/dP slope at room temperature. 

Since the sc phase of P becomes superconducting (SC) 
below 10 K, one may expect similar electron-phonon 
interaction (EPI) strengths in P at higher pressures. 
Hence, a conventional s-wave pairing state may also arise 
in bcc P. We investigate this possibility by using a linear- 
response full-potential LMTO DFT method |2j, which 
can describe the SC properties for sp materials with 
high accuracy |22l as well as those of d-metals under 
high pressure. [23| The calculated FP-LMTO electronic 
structure of bcc P is practically identical to that of the 



WIEN-code calculation. The phonon dispersions w m (q) 
and density of states F(u>) are obtained from this linear- 
response technique [2l| and used to calculate the Eliash- 
berg spectral function a 2 (oj)F(oj) which defines the EPI 



strength A = J * 



OJ 



- a ■ (u>) F(<jj) duj. The function 



a 2 (uj)F(ui), shown in Fig. 2, is dominated by the two 
optical peaks at 2.1 THz and 8.3 THz. 

The SC transition temperature is calculated by means 
of the Allen-Dynes modified McMillan expression j2(J 

using the calculated EPI 



T c = 



1.2 



exp 



-1.04(1+A) 



A-^*(1+0.62A) 

strength and assuming a value of the Coulomb pseudopo- 
tential of /x*=0.15 which was also used to calculate T c for 
the sc phase of P. In bcc P we find T c =19 K. Regarding 
the transport properties of bcc P in the normal state, 
the temperature dependence of the specific dc resistivity 
shows a crossover from power law to linear dependence 
at 50 K. 

There is considerable current interest in technological 
spintronics applications, in which the spin of carriers is 
exploited to provide new functionality for microelectronic 
devices. Growth of P containing materials by thin film 
techniques, such as molecular beam epitaxy or pulsed 
laser deposition, has been developed to fabricate the class 
of wide gap phosphide ferromagnetic semiconductors . [Toj 
Epitaxial thin film structures offer ideal opportunities to 
stabilize the bcc phase of P at ambient conditions using 
suitable growth templates such as the V(100), Fe(100) or 
Cr(100) substrates. Consequently, we now turn to study 
the electronic and magnetic properties of thin films of 
P (3-5 layers in thickness) embedded in ferromagnetic 
(FM) bcc-Fe. During the last decade considerable theo- 
retical and experimental progress has been made towards 
the understanding of metallic magnetism in superlattices 
and multilayers. However, to date, no such investigations 
have been made for Fe/P. 

The electronic structure of the Fe5/P„/Fe4 (n=3,5) 
superlattices, with the lattice constant of bulk Fe, is cal- 
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FIG. 3: Total and layer-resolved electronic DOS of the 
Fe5/Pa/Fe4 superlattices with the lattice constant of bulk Fe. 
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FIG. 4: Calculated layer-resolved spin moments of Fe and 
P in the Fe5/P„/Fe4 (n=3,5) superlattices with the lattice 
constant of bulk Fe. 



culated from the WIEN-code. Fig. 3 shows the total 
and layer-resolved electronic density of states (DOS) for 
n = 3. Evidently the DOS of phosphorus is fairly insen- 
sitive to monolayer (ML) position and p-d hybridization 
has broadened d-related peaks in the Fe interface DOS. 
Following on from this Fig. 4 shows the layer-resolved 
spin magnetic moments and the spin magnetic moment 
per atom on the Fe interface layer is about 20 % smaller 



than that of bulk bcc Fe whereas the second Fe layer 
from the ideally sharp Fe/P interface has an enhanced 
local spin moment compared to that of bulk Fe. The P 
atoms away from the interface have no induced magnetic 
moments. The spin magnetic moment on the P interface 
atoms are tiny and antiparallel to the Fe magnetization 
direction. Thus, the total magnetic moments of the 
Fe/P/Fe superlattices, which are almost independent of 
the thickness of P, are slightly reduced from those of a 
relevant pure Fe system due to the AF Fe-P interface 
coupling. Obviously, this effect would be important for 
superlattices with thin Fe slabs. We note that a simi- 
lar but more pronounced AF coupling takes place at the 
Fe/V (001) interface. Hill 

The proximity effect in FM/SC heterostructures means 
that the Cooper pair amplitude in an exchange field does 
not decay exponentially to zero on the FM side but oscil- 
lates with decreasing amplitude as a function of the dis- 
tance from the interface. These oscillations occurring in 
what is known as the Fulde-Ferrell-Larkin-Ovchinnikov 
state m may lead striking changes in the SC transi- 
tion temperature (as observed in Fe/V and Fe/Nb het- 
erostructures |23) as a function of the FM slab thickness 
d,Fe- We suggest that Fe/P heterostructures may show 
similar effects. 

Another interesting application of the SC state in the 
stabilized bcc phase of P might be available with V/Fe/P 
junctions, in which thin FM films are sandwiched be- 
tween two superconductors with differing SC transition 
temperatures. Such systems support spontaneous cur- 
rents parallel to the FM/SC interface while the spin po- 
larization of these currents depends on band filling [29| 
and, hence, can be readily adjusted. 

In summary, our study of P under high pressure based 
upon DFT and Debye model calculations produces a 
pressure-temperature phase diagram in good agreement 
with experiment. Pressure-induced metallization, at- 
tributed to the strong P-P interactions, is enhanced in 
the bcc phase which is stable for pressures in excess of 
265 GPa. We have demonstrated that bcc-P might be 
a s-wave superconductor below 19 K owing primarily to 
strong EPI for the phonon bands along [110]. Finally, we 
predict that the SC properties of bcc-P can be realized in 
Fe/P heterostructures in which the bcc phase of P should 
be stable. 
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